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ABSTRACT 

Aims. The Herschel-ATLAS survey (H-ATLAS) will be the largest area survey to be undertaken by the Herschel Space Observatory. It will cover 

550 sq. deg. of extragalactic sky at wavelengths of 100, 160, 250, 350 and 500/jm when completed, reaching flux limits (5cr) from 32 to 145mJy. 

We here present galaxy number counts obtained for SPIRE observations of the first ~14 sq. deg. observed at 250, 350 and 500//m. 

Methods. Number counts are a fundamental tool in constraining models of galaxy evolution. We use source catalogs extracted from the H-ATLAS 

maps as the basis for such an analysis. Correction factors for completeness and flux boosting are derived by applying our extraction method to 

model catalogs and then applied to the raw observational counts. 

Results. We find a steep rise in the number counts at flux levels of 100 — 200mJy in all three SPIRE bands, consistent with results from BLAST. 

The counts are compared to a range of galaxy evolution models. None of the current models is an ideal fit to the data but all ascribe the steep rise 

to a population of luminous, rapidly evolving dusty galaxies at moderate to high redshift. 

Key words. Galaxies: evolution; Galaxies: statistics; Infrared: galaxies; Submillimetre: galaxies 
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1. Introduction 

The determination of the number of sources observed as a func- 
tion of brightness is fundamental in observational astronomy. 
These so-called number counts have been used in the past to 
determine the structure of our Galaxy (Herschel, 1784), rule 
. out cosmological models (Ryle & Clark, 1961), and examine 

■ the evolution of galaxies and quasars (eg. Small et al., 1995; 
' Hasinger et al., 1993). Number counts have revealed that very 

rapid evolution must take place among infrared luminous popu- 
lations to account for the large number of sources seen at faint 
, fluxes in the submillimetre (eg. Hughes et al., 1998; Small et 

■ al., 1997; Bales et al., 2000) and to produce the cosmic infrared 
I background (Puget et al., 1996; Fixsen et al., 1998; Hauser & 

Dwek 2001). The detailed origin and nature of this evolution re- 
mains unclear following observations in the mid-IR (eg. Elbaz et 
al., 1999, Papovich et al., 2004, Le Floc'h et al., 2009; Babbedge 
1 et al., 2006) and far-IR (eg. Frayer et al., 2006) using ISO and 
Spitzer, and extensive followup studies of submm sources (eg. 
Chapman et al., 2005; Clements et al., 2008; Dye et al., 2008). 
Issues needing to be addressed include spectral energy distribu- 
tion (SED) evolution and the redshift distribution of the evolving 
luminous populations. Understanding these processes is a key 
goal of the Herschel Space Observatory (Pilbratt et al., 2010) 
and of the various large extragalactic surveys that it is under- 
taking. We here report the first number counts obtained by the 
Herschel-ATLAS (H-ATLAS) survey (Bales et al., 2010) which 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



set new constraints on models of galaxy evolution in the far- 
IR/submm. 



2. Observations and Source Extraction 

The overall H-ATLAS survey is described by Bales et al. (2010). 
Briefly, it is a survey that will cover ~550 sq. deg. of extra- 
galactic sky at wavelengths of 100, 160, 250, 350 and 500/vm 
to 5cr depths of 32-145mJy. The observations discussed here 
cover a ~14 sq. deg. region, <3% of the total area, which 
were obtained during the Science Demonstration Phase (SDP) 
of Herschel operations. This field is centred on RA=09:05:30.0 
DBC=00:30:00.0. The field lies in the centre of one of the 
Galaxy And Mass Assembly (GAMA) optical survey fields 
(Driver et al., 2009), GAMA9, so plentiful complementary data 
from UV to near-IR, including optical spectra, is available. 

H-ATLAS uses parallel mode observations which provide 
data simultaneously from both the PACS (Poglitsch et al. 2010) 
and SPIRB (Griffin et al. 2010) instruments. The time-line data is 
reduced using a naive mapping technique after removing instru- 
mental temperature variations from the time-line data (Pascale 
et al. in prep., see also Griffin et al., 2008). Noise maps were 
generated using the two cross-scan measurements to estimate the 
noise per detector pass, and then for each pixel the noise is scaled 
by sqrt(number of detector passes). PACS data and results will 
be discussed elsewhere. 

Sources were identified in the SPIRB maps by first subtract- 
ing a local background, estimated from the peak of the histogram 
of pixel values in 30 x 30 blocks, corresponding to 2.5' for the 
250jum map, and 5' for the 350 and 500jum maps. Larger blocks. 
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with 60 X 60 and 120 x 120 pixels were tried, but these pro- 
duced no significant differences (>~ Icr) in derived counts. The 
background at each pixel was then estimated using a bi-cubic 
interpolation between the coarse grid of backgrounds. The back- 
ground subtracted maps were then filtered by the estimated PSF, 
including inverse variance weighting, where the noise for each 
pixel was estimated from the noise map. 

The maps from all three bands are then combined with 
weights set by the local inverse variance, and also the prior ex- 
pectation of the SED of the galaxies. We tried a flat-spectrum 
prior, where equal weight is given to each band and also 250yum 
weighting, where only the 250//m was included. At the depth of 
the filtered maps source confusion is a significant problem, and 
the higher resolution of the 250yum maps outweighed the signal- 
to-noise gain from adding in the other bands. So for our current 
data we chose to use the 250/im only prior for all our catalogues. 
We plan to revisit this issue in future data-releases, but also com- 
pared the number counts produced by our chosen method and the 
flat spectrum prior and found no significant differences. Further 
information on our background subtraction and source extrac- 
tion methods can be found in Maddox et al. (2010) and Rigby et 
al. (in preparation). 

All local peaks are identified in the 250yum PSF filtered map 
as potential sources, and a Gaussian is fitted to each peak. This 
provides an estimate of the position at the sub-pixel level, and an 
estimate of the peak value, which gives the best flux estimate for 
a point source. The fluxes in other bands were estimated by using 
a bi-cubic interpolation to the position given by the 250/im map. 
This extraction method was validated against a range of source 
fluxes in our simulations, and will be discussed in detail in Rigby 
et al. (in prep). To produce a catalogue of reliable sources, we se- 
lected only sources that are detected at the 5-cr level in any of the 
bands. In calculating the cr for each source, we use the relevant 
noise map, and add the confusion noise to this in quadrature. 
The average 1-cr instrumental noise values for the PSF filtered 
maps are 4, 4 and 5.7mJy/beam respectively in the 250, 350 and 
500/im bands. We estimated the confusion noise from the differ- 
ence between the variance of the maps and the expected variance 
due to instrumental noise, and find that the 1-cr confusion noise 
is 5, 6 and 7 mJy/beam at 250, 350 and 500/im. The resulting to- 
tal 5-cr limits are 33 36 and 45mJy/beam. It should also be noted 
that overall flux calibration uncertainties are at the 15% level for 
SPIRE (Swinyard et al., in prep). 
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Table 1. H- ATLAS Galaxy Counts. Correction column gives the 
effective flux correction factor required to match the output to 
input counts from the simulations, and the fluxes quoted are after 
correction. 



3. The H-ATLAS Number Counts 

We use our extracted catalogs as the basis for the determi- 
nation of the H-ATLAS number counts in the SPIRE bands. 
Additionally, those objects identified with optical galaxies 
whose size extends beyond the instrumental beam, have fluxes 
from an explicit aperture extraction with size equal to the in- 
strument beam added in quadrature to the measured object size. 
One further source, a merger of two galaxies where the opti- 
cal radius is too small to include the second component, had an 
aperture assigned by hand. An additional very bright, extended 
source, known to be galactic in origin (Thompson et al., in prep), 
was excluded from our analysis. The resulting catalog of source 
fluxes in the three bands is cut so that only sources detected at 
5cr or above in a given band are used for the counts. 

Completeness, reliability and the effects of both flux boost- 
ing (eg. Coppin et al. 2006) and blending were assessed using 
simulations of our observations. These have the same noise prop- 
erties as the processed maps and include a realistic cirrus back- 
ground, based on IRAS measurements. Simulated sources were 



generated using the model of Negrello et al (2007), with an ad- 
ditional reduction in flux density of 15% to provide better agree- 
ment with the measured counts. The source extraction method 
described above was applied to these images, and the resulting 
simulated catalogue was also cut at 5cr in each band. Flux cor- 
rection factors matching input and output counts, and account- 
ing for the effects of both boosting/blending and incompleteness, 
were derived from the simulations. Where greater than 5%, these 
corrections were applied to the real counts to produce our final 
numbers. The band mainly affected by this is the 500/im channel 
where flux boosting and blending are most significant due to the 
large beams. There is also a correction to the lowest flux bin of 
the 250//m band as a result of incompleteness. 

The Euclidean normalised number counts derived from this 
data and using the corrections shown in Table 1 are shown in 
Fig. 1, alongside a variety of models and data from BLAST 
(Patanchon et al., 2009) and, at 350/im, Khan et al. (2007). Table 
1 also provides raw counts, normalised differential counts and 
integral counts. 
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4. Discussion 





10 


m 








>> 

-3 


10 


CO 




O) 








, — , 




\j 


10 


o 




lO 




C\J 




CO 








■o 


10 



" 10^ 






o 

in 

CO 
CO 

W 
"a 



0.01 



0.01 



0.01 



. 


"^ 


J*- 




= 


H-ATLAS 

BLAST 

Lagache 


• 




- 


Lacey 
Willman — 


: 










:^ 


=^ 









0.10 
F250 (Jy) 



1.00 



in' 


; \ 











' ' ' ': 








H-ATLAS 

BLAST 

SHARC2 

Lagache 

Pearson — 


• 


10^ 


^5^ 


^ i ^^^ 




: 




Willman — 


— : 






_ 


10^ 




%. 

X 


^ 


f^ 












10^ 









0.10 
F350 (Jy) 



1.00 




0.10 
F500 (Jy) 



1.00 



Fig. 1. H-ATLAS Euclidean Normalised Differential Number 
Counts compared to models and other data. Data from this work 
are shown as solid dots. Data from BLAST is shown as open 
diamonds. The solid horizontal line shows the asymptotic lo- 
cal differential counts from Serjeant & Harrison (2005) with 
dashed lines showing the ±11% fractional uncertainly. Models 
are shown as follows: Lagache et al. (2004) - green; Pearson 
& Khan (2009) - red; Negrello et al. (2007) - blue; Lacey et 
al.(2008) - orange; Wilman et al. (2010) - purple. Top: 250yum 
counts. Middle: 350;um counts, with the 350//m point of Khan 
et al. (2007) shown as an open triangle. Bottom: SOOjum counts. 



We compare the H-ATLAS counts to a range of models and to 
existing data from BLAST (Patanchon et al., 2009), SHARC2 
(Khan et al., 2007) at 350yum, and the Serjeant & Harrison (2005) 
study of the local far-IR/submm luminosity function. The mod- 
els include those of Lagache et al (2004), Pearson et al. (2009), 
Negrello et al. (2007), Lacey et al., (2009) and Wilman et al. 
(2010) and cover a wide range of modelling techniques and 
assumptions. The Lagache et al. (2004) model is a backward- 
evolution model in which two template SEDs, one for nor- 
mal galaxies and another for starbursts, are evolved indepen- 
dently from z=0. The Pearson et al. (2009) model is a simi- 
lar backward-evolution model but instead uses seven luminosity 
dependent template SEDs. The rate of evolution is also lumi- 
nosity dependent with stronger evolution in the more luminous 
sources (bright mode evolution in the terminology of Pearson 
et al. (2009)). The Negrello et al. (2007) count model includes 
contributions from a range of populations. These include submm 
galaxies evolving according to a physical model (Granato et al., 
2004), radio sources, as modelled by DeZotti et al. (2005), and 
late-type (starburst plus normal spiral) galaxies. Unlike the oth- 
ers, the counts of SMGs in this model also include a contribution 
from strongly-lensed galaxies, quantified following the recipes 
of Perrotta et al. (2002). The Lacey model is calculated from a 
theoretical model of galaxy formation based on CDM (Lacey et 
al. 2010). It uses the GALFORM semi-analytical galaxy forma- 
tion model in combination with the GRASIL model for the re- 
processing of starlight by dust (Granato et al 2000). The galaxy 
formation model is the same as that described in Baugh et al 
(2005), and compared to Spitzer data in Lacey et al (2008), and 
incorporates a top-heavy IMF in starbursts triggered by galaxy 
mergers. The Wilman et al. (2010) source counts are based on 
the large SKADS radio continuum simulation (Wilman et al. 
2008) and uses the radio-FIR relation together with a model of 
the evolution of this relation, along with various prescriptions 
to account for the AGN and starburst populations, to predict the 
counts in the FIR/submm. All of these models are broadly con- 
sistent with existing data on galaxy counts, redshift distributions 
and luminosity functions from the mid-IR to radio bands, includ- 
ing data from IRAS, Spitzer and ground-based submm surveys. 

A clear feature in the observed counts is a steep increase 
above an extrapolation of the bright counts to fainter fluxes. This 
increase starts at fluxes -100 — 200mJy. Despite being hinted 
at in earlier observations, Herschel allows this rise in counts 
to be unambiguously detected and studied in detail for the first 
time. The feature is seen in all three SPIRE bands, but is most 
pronounced in the longest wavelength, 350 and 500/im, bands 
where the counts rise roughly twice as much as at 250/^m. Such 
a feature is predicted in a number of models though no model is 
a particularly good fit to the counts in all bands. The best per- 
forming model with respect to the shape of the counts is that 
of Negrello et al. (2007), but this model overpredicts the effect 
of the bump at 250 and 350 yum. The Negrello model also pro- 
vides a better fit to the brighter counts than most of the other 
models, especially at 350/im, but there is clearly much work 
to be done to improve the match between any of these mod- 
els and the observed data. The number count rise represents the 
strongly evolving luminous (> lO^L©) galaxy population that is 
a generic feature in all galaxy evolution models that are able to fit 
the submm galaxies. The details of this evolution, possibly com- 
bined with other effects such as gravitational lensing (a feature of 
the Negrello model), control the shape and position of the rise. 
Its strength as a function of wavelength will be dependent on 



D.L. Clements et al.: The Herschel- ATLAS: Extragalactic Number Counts from 250 to 500 Microns 



the detailed redshift and possibly evolving SEDs of the rapidly- 
evolving population which can also be constrained by the colour 
of individual sources (see eg. Amblard et al. 2010). Its promi- 
nence at 350 and SOOjum compared to 250jum thus argues that 
the sources responsible lie at z>1.5, as suggested by Granato 
et al. (2004). If they lie at lower redshifts they would require 
colder SEDs. It is also interesting to note that some models (eg. 
that of Lacey) appear to be strongly excluded by the observations 
and that the largest discrepancy between predicted and observed 
counts is at the bright end, where the full H- ATLAS survey will 
be uniquely powerful. These models are also unable to repro- 
duce the asymptotic local counts from Serjeant and Harrison 
(2005) based on IRAS data. This suggests that our knowledge 
of brighter, more local far-IR/submm galaxies is not as good as 
we might wish. 

These new counts from H- ATLAS set the scene for a much 
better understanding of the enigmatic far-IR/submm population. 
The SDP data presented here represents less than 3% of the total 
area that will be covered by H- ATLAS. Although this is only 
a small fraction of H-ATLAS, these data already include more 
than 5 times as many (> 5cr) sources as previously published 
at these wavelengths. As further data is added, more detailed 
measurements towards the bright end will allow the effects of 
strong lensing and local sources on the counts to be determined. 
H-ATLAS is also large enough to totally overcome the effects of 
cosmic variance, and allow it to be studied by comparing counts 
in sub-regions. 

5. Conclusions 

We have used the first ~14 sq. deg. of the H-ATLAS project 
observed by Herschel to determine the number counts of far- 
IR/submm sources in the SPIRE bands at 250, 350 and 500/vm. 
We find that there is a steep rise in the counts at fluxes of 100- 
200mJy. These observations are compared to a range of galaxy 
evolution models and data from BLAST and elsewhere. The H- 
ATLAS counts are found to be consistent with previous observa- 
tions, but no theoretical model is an ideal fit. Some of the models 
are quite poor matches to the observations, with counts at fluxes 
brighter than ~200mJy setting strongly excluding some models. 
The best fitting model of those considered is that of Negrello 
et al. (2007). These observations represent just 3% of the total 
area to be covered by the H-ATLAS survey. Once completed this 
project will clearly be a very powerful tool for understanding the 
far-IR/submm galaxy population. 
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